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Abstract. Lunar Laser Ranging (LLR) measurements are crucial for advanced exploration of 
the evolutionary history of the lunar orbit, the laws of fundamental gravitational physics, seleno- 
physics and geophysics as well as for future human missions to the Moon. Current LLR technique 
measures distance to the Moon with a precision approaching 1 millimeter that strongly demands 
further significant improvement of the theoretical model of the orbital and rotational dynamics 
of the Earth-Moon system. This model should inevitably be based on the theory of general rela- 
tivity, fully incorporate the relevant geophysical/ selenophysical processes and rely upon the most 
recent I A U standards. We discuss new methods and approaches in developing such a mathematical 
model. The model takes into account all classic and relativistic effects in the orbital and rotational 
motion of the Moon and Earth at the millimeter-range level. It utilizes the lA U 2000 resolutions on 
reference frames and demonstrates how to eliminate from the data analysis all spurious, coordinate- 
I dependent relativistic effects playing no role in selenophysics/ geophysics. The new model is based 

on both the locally-inertial and barycentric coordinates and extends the currently used LLR code. 
The new theory and the millimeter LLR will give us the opportunity to perform the most precise 
O , fundamental test of general relativity in the solar system in robust and physically- adequate way. 

Lunar Laser Ranging 

The dynamical modeling for the solar system (major and minor planets), for deep space navi- 
[ gation, and for the dynamics of Earth's satellites and the Moon must be consistent with general 

I relativity. Lunar laser ranging (LLR) measurements arc particularly crucial for testing general rcl- 

~l ' ativistic predictions and advanced exploration of other laws of fundamental gravitational physics. 

, Current LLR technologies allow us to arrange the measurement of the distance from a laser on the 

Earth to a corner-cube reflector (CCR) on the Moon with a precision approaching 1 millimeter 
[1, 2]. There is a proposal to place a new CCR array on the Moon [3], and possibly to install 
other devices such as microwave transponders [4] for multiple scientific and technical purposes. 
I Successful human exploration of the Moon strongly demands further significant improvement of 

the theoretical model of the orbital and rotational dynamics of the Earth-Moon system. This 
model should inevitably be based on the theory of general relativity, fully incorporate the relevant 
. geophysical processes, lunar libration, tides, and should rely upon the most recent standards and 

^ ' recommendations of the lAU for data analysis [5] . 

5^ , LLR technique is currently the most effective way to study the interior of the Moon and dynamics 

of the Moon-Earth system. The most important contributions from LLR include: detection of a 
molten lunar core and measurement of its influence on the Moon's orientation along with tidal 
dissipation [6-9]; detection of lunar free libration along with the forced terms from Venus [10] and 
the internal excitation mechanisms [11]; an accurate test of the strong principle of equivalence 
for massive bodies [12, 13] also known as the Nordtvedt effect [14, Section 8.1]; and setting of a 
stringent limit on time variability of the universal gravitational constant and (non)cxistencc of long- 
range fields besides the metric tensor [15]. LLR analysis has also given access to more subtle tests 
of relativity [16-19], measurements of the Moons tidal acceleration [20-22] and geodetic precession 
of the lunar orbit [23, 24], and has provided orders-of- magnitude improvements in the accuracy of 
the lunar ephemeris [25-29] and its three-dimensional rotation [30, 31]. On the geodesy front, LLR 
contributes to the determination of Earth orientation parameters, such as nutation, precession 
(including relativistic geodetic precession), polar motion, UTl. and to the long-term variation of 
these effects [32, 33]. 

LLR also contributes to the realization of both the terrestrial and sclcnodesic reference frames [34, 
35]. The satellite laser ranging (SLR) realization of a dynamically-defined inertial reference frame 
[36] in contrast to the kinematically-realized frame of VLBI [37, Section 6] , offers new possibilities 
for mutual cross-checking and confirmation [33] especially after the International Laser Ranging 
Service (ILRS) was established in September 1998 to support programs in geodetic, geophysical, 
and lunar research activities and to provide the International Earth Rotation Service (lERS) with 
products important to the maintenance of an accurate International Terrestrial Reference Frame 
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(ITRF) [38]. 

Over the years, LLR has benefited from a number of improvements both in observing technol- 
ogy and data modehng [39]. Recently, sub-centimeter precision in determining range distances 
between a laser on Earth and a retro-reflector on the Moon has been achieved [1, 2]. As precision 
of LLR measurements was gradually improving over years, enormous progress in understanding 
evolutionary history of the Earth-Moon orbit and the internal structure of both planets has been 
achieved. With the precision approaching 1 millimeter and better, accumulation of more accurate 
LLR data will lead to new, fascinating discoveries in fundamental gravitational theory, geophysics, 
and physics of lunar interior [40] whose unique interpretation will intimately rely upon our ability 
to develop a systematic theoretical approach to analyze the sub-centimeter LLR data [28, 41]. 



EIH Equations of Motion in N-body Problem 

Nowadays, the theory of the lunar motion should incorporate not only the numerous Newto- 
nian perturbations but has to deal with much more subtle relativistic phenomena being currently 
incorporated to the ephemeris codes [25, 42-44]. Theoretical approach, used for construction of 
the ephemerides, accepts that the post-Newtonian description of the planetary motions can be 
achieved with the Einstein-Infeld-Hoffmann (EIH) equations of motion of point-like masses [45], 
which are valid in the barycentric frame of the solar system with time coordinate, t, and spatial 
coordinates, = x. 

Due to the covariant nature of general theory of relativity the barycentric coordinates are not 
unique and are defined up to the space-time transformation [46-48] 
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where summation goes over all the massive bodies of the solar system {B = 1,2, A^); G is the 
universal gravitational constant; c is the fundamental speed in the Minkowskian space-time; a dot 
between any spatial vectors, a ■ b, denotes an Euclidean dot product of two vectors a and b; M^, 
is mass of a body B; Xb ~ Xsit) and Vb = Vsit) are coordinates and velocity of the center of 
mass of the body B; Rg relative distance from a field point x to the body B; i/b 

and Xb are constant, but otherwise free parameters being responsible for a particular choice of the 
barycentric coordinates. We emphasize that these parameters can be chosen arbitrary for each 
body B of the solar system. Physically, it means that the space-time around each body is covered 
locally by its own coordinate grid, which matches smoothly with the other coordinate charts of 
the massive bodies in the buffer domain, where the different coordinates overlap. 

If the bodies in N-body problem are numbered by indices B,C,D, etc., and the coordinate 
freedom is described by equations (l)-(2), EIH equations of motion for the body B have the 
following form [46] 
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where E]^^ is called the gravitoelectric force, and the terms associated with the cross products 
{vb X HgcT and (vc x HbcY are referred to as the gravitomagnetic force [49]. 
The gravitoelectric force is given by 



Ebc - ( — ST^ ^ T~^Bc 1 Rbc 1 (4) 
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and the gravitomagnetic force is 
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where Vbc = Vb —Vq. The first term in right side of equation (4) is the Newtonian force of gravity. 



3 



and the post-Newtonian correction 

£bc = -^^{~l + \c)vl+^{l-2\c){vs-Vc)-{l-Z\c)vl (6) 
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As one can see, the gravitomagnctic force (5) is proportional to the gravitoelectric force (4) mul- 
tiplied by the factor of V/c, where V is the relative velocity between two gravitating bodies. EIH 
equations (3)-(5) differ from the equations of the PPN formalism 1969ApJ...158...81E employed in 
particular at JPL for actual calculation of the ephemerides of the major planets by the fact that 
the right side of equation (3) has been resolved into radius-vectors and velocities of the massive 
bodies and does not contain second derivatives (accelerations). This elimination of the high-order 
time derivatives from a perturbed force is a standard practice in celestial mechanics for calculation 
of the perturbed motion. 

Barycentric coordinates Xb and velocities Vb of the center of mass of body B are adequate 
theoretical quantities for description of the world-line of the body with respect to the center of 
mass of the solar system. However, the barycentric coordinates are global coordinates covering the 
entire solar system. Therefore, they have little help for efficient physical decoupling of the post- 
Newtonian effects existing in the orbital and rotational motions of a planet and for the description 
of motion of planetary satellites around the planet. The problem stems from the covariant nature 
of EIH equations, which originates from the fundamental structure of space-time manifold and the 
gauge freedom of the general relativity theory. 

This freedom is already seen in the post-Newtonian EIH force (3)-(6) as it explicitly depends on 
the choice of spatial coordinates through parameters Ac , Ad . Each term depending explicitly on 
Ac, Ad in equations (3)-(6), has no direct physical meaning as it can be eliminated after making a 
specific choice of these parameters. In many works on experimental gravity and applied relativity 
researches fix parameters Ac = Ad = 0, which corresponds to working in harmonic coordinates. 
Harmonic coordinates simplify EIH equations to large extent but one has to keep in mind that 
they have no physical privilege anyway, and that a separate term or a limited number of terms 
from EIH equations of motion can not be measured - only coordinate-independent effects can be 
measured [46]. 

Recently, there was a lot of discussions about whether LLR can measure the gravitomagnctic field 
-ff^c [1^7 19' 51-53]. The answer to this question is subtle and requires more profound theoretical 
consideration involving the process of propagation of the laser pulses in a curved space-time of the 
Earth-Moon system. We are hoping to discuss this topic somewhere else. Nevertheless, what is 
evident already now is that equation (3) demonstrates a strong dependence of the gravitomagnctic 
force on the choice of coordinates. For this reason, by changing the coordinate parameter Ac one 
can eliminate either the term {vb x HbcY or {vc x HbcY from EIH equations of motion (3). It 
shows that the strength of the factual gravitomagnctic force, as it appears in the equations of 
motion, is coordinate-dependent, and, hence, a great care should be taken in order to properly 
interpret the LLR " measurement" of such gravitomagnctic terms in consistency with the covariant 
nature of the general theory of relativity and the theory of astronomical measurements in curved 
space-time [48, 54-57]. 



The Gauge Freedom 



The primary gauge freedom of EIH equations of motion is associated with the transformations 
(l)-(2) of the barycentric coordinates of the solar system, which are parameterized by parameters 
i>c a-nd Ac. However, the post-Newtonian force in the lunar equations of motion admits addi- 
tional freedom of coordinate transformations. This residual freedom remains even after fixing the 
coordinate parameters Vc and Ac in equations (3)-(6). It is associated with the fact that the 
Earth-Moon system moves in tidal gravitational field of the Sun and other planets, which presence 
in the local frame of the Earth-Moon system indicates that the local background space-time is 
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not asymptotically-flat. The residual freedom remains in making transformations of the local co- 
ordinates attached to the Earth-Moon system. It induces the gauge transformation of the metric 
tensor and the Christoffel symbols and changes the structure of the post-Newtonian terms in EIH 
equations of motion of the Earth-Moon system that is not associated with the parameters Vc Stiid 
Ac. 

Thus, we face the problem of investigation of the residual gauge freedom of the lunar equations 
of motion, which goes beyond the choice of the barycentric coordinates by fixing a specific value 
of the gauge parameter Ac in equations (3)-(6). This freedom is naturally associated with the 
choice of the local coordinates of the Earth-Moon barycentric frame as well as the geocentric and 
selenocentric reference frames. Proper choice of the local coordinates removes all non-physical 
degrees of freedom from the metric tensor and eliminates spurious (non-measurable) terms from 
the post-Newtonian forces in the equations of relative motion of the Moon with respect to the 
Earth. If one ignores the residual gauge freedom, the gauge-dependent terms will infiltrate the 
equations of motion causing possible misinterpretation of LLR observations. This problem is 
similar to that one meets in cosmology, where the theory of cosmological perturbations is designed 
essentially in terms of the gauge-independent variables so that observations of various cosmological 
effects are not corrupted by the spurious, coordinate-dependent signals [58]. The residual gauge 
degrees of freedom existing in the relativistic three-body problem (Sun-Earth-Moon), can lead to 
misinterpretation of various aspects of gravitational physics of the Earth-Moon system [51, 59], 
thus, degrading the value of extremely accurate LLR measurements for testing fundamental physics 
of space-time and deeper exploration of the lunar interior [28] . 

The residual gauge freedom of the three body problem was studied by Brumberg and Kopeikin 
[60], Klioner and Voinov [61], and Damour, Soffel and Xu [62]. They found that the post-Newtonian 
equations of motion of a test body (artificial satellite) can be significantly simplified by making use 
of a four-dimensional space-time transformation from the solar barycentric coordinates = {ct, x), 
to the geocentric coordinates X" = {cT,X) 



where the gauge functions A{t,x), B{t,x), C'^{t,x) are polynomials of the geocentric distance 
Te = X ~ XE{t) of the field point x from the Earth's geocenter, XE{t). Coefficients of these 
polynomials are functions of the barycentric time t that are determined by solving a system of 
ordinary differential equations, which follow from the gravity field equations and the tensor law of 
transformation of the metric tensor from one coordinate chart to another [63] . Contrary to the test 
particle, the Moon is a massive body, which makes the exploration of the residual gauge freedom 
of the lunar motion more involved. This requires introduction of one global (SSB) frame and three 
local reference frames associated with the Earth-Moon barycenter, the geocenter, and the center of 
mass of the Moon (selenocenter). It should be clearly understood that any coordinate system can 
be used for processing and interpretation of LLR data since any viable theory of gravity obeys the 
Einstein principle of relativity, according to which there is no preferred frame of reference [64-66] . 
Accepting the Einstein principle of relativity leads to discarding any theory of gravity based on a 
privileged frame (aether) [67] or admitting a violation of the Lorentz invariance [68]. The class of 
scalar-tensor theories of gravity, which have two PPN parameters - /3 and 7 [14, 69], is in agreement 
with the principle of relativity and it has been investigated fairly well [70, 71]. 

The principle of relativity also assumes that an arbitrary chosen, separate term in the post- 
Newtonian equations of motion of massive bodies can not be physically interpreted as straightfor- 
ward as in the Newtonian physics. The reason is that the post-Newtonian transformations (l)-(2) 
and (7)-(8) of the barycentric and local coordinates, change the form of the equations of motion 
so that they are not form- invariant. Therefore, only those post-Newtonian effects, which do not 
depend on the frame transformations can have direct physical interpretation. For example, the 
gauge parameters Vc and Ac entering transformations (l)-(2) and EIH equations (3)-(5) can not 
be determined from LLR data irrespectively of their accuracy because these parameters define the 
barycentric coordinates and can be fixed arbitrary by observer without any relation to observations. 




(8) 



Towards a New Lunar Ephemeris 



Existing computer-based theories of the lunar ephemeris [25, 42-44] consist of three major blocks: 



5 



(1) the baryccntric EIH equations (3)-(6) of orbital motion of the Moon, Earth, Sun, and other 
planets of the solar system with the gauge parameters Vc ~ 1/2, Ac = 0- the standard PPN 
coordinates; 

(2) the Newtonian rotational equations of motion of the Moon and Earth; 

(3) the barycentric post-Newtonian equations of motion for light rays propagating from laser to 
CCR on the Moon and back in standard coordinates with the gauge parameters Vc = 1/2, 
Ac = 0. 

This approach is straightforward but it does not control gauge-dependent terms in EIH equations 
of motion associated with the choice of the gauge-fixing parameters Uc and Ac • Particular disad- 
vantage of the barycentric approach in application to the lunar ephemeris is that it mixes up the 
post-Newtonian effects associated with the orbital motion of the Earth-Moon barycenter around 
the Sun with those, which are attributed exclusively to the relative motion of the Moon around 
the Earth. This difficulty is also accredited to the gauge freedom of the equations of motion in 
three-body problem and was pointed out in papers [60, 62, 72]. Unambiguous decoupling of the 
orbital motion of the Earth-Moon barycenter from the relative motion of the Moon around the 
Earth with apparent identification of the gauge-dependent degrees of freedom in the metric ten- 
sor and equations of motion is highly desirable in order to make the theory more sensible and to 
clean up the LLR data processing software from the fictitious coordinate-dependent perturbations, 
which do not carry out any physically-relevant information and may accumulate errors in numerical 
ephemerides of the Moon and the Earth. 

This goal can be rationally achieved if the post-Newtonian theory of the lunar motion is con- 
sistently extended to account for mathematical properties offered by the scalar-tensor theory of 
gravity and the differential structure of the space-time manifold. Altogether it leads us to the 
idea that besides the global barycentric coordinates of the solar system one has to introduce three 
other local reference frames. The origin of these frames should be fixed at the Earth-Moon system 
barycenter, the Earth's center of mass (geocenter), and the Moon's center of mass (selenocenter) . 
We distinguish the Earth-Moon barycenter from the geocenter because the Moon is not a test 
particle, thus, making the Earth-Moon barycenter displaced from the geocenter along the line con- 
necting the Earth and Moon and located approximately 1710 km below the surface of the Earth. 
Mathematical construction of each frame is reduced to finding a metric tensor by means of solution 
of the gravity field equations with an appropriate boundary condition [41]. The gauge freedom of 
the three-body problem is explored by means of matching the set of the metric tensors defined in 
each reference frame in the overlapping domains of their applicability associated with the specific 
choice of boundary conditions imposed in each frame on the metric tensor. This matching proce- 
dure is an integral part of the equations defining the local differential structure of the manifold 
[73, 74], which proceeds from a requirement that the overlapping space-time domains covered by 
the local reference frames, are diffeomorphic. 

The primary objective of the multi-frame post-Newtonian theory of the lunar ephemeris is the 
development of a new set of analytic equations to revamp the LLR data processing software in order 
to suppress the spurious gauge-dependent solutions, which may overwhelm the existing barycentric 
code at the millimeter accuracy of LLR measurements, thus, plunging errors in the interpretation of 
selenophysics, geophysics and fundamental gravitational physics. Careful mathematical construc- 
tion of the local frames with the post-Newtonian accuracy will allow us to pin down and correctly 
interpret all physical effects having classical (lunar interior. Earth geophysics, tides, asteroids, etc.) 
and relativistic nature. The gauge freedom in the three-body problem (Earth-Moon-Sun) has been 
carefully examined in our paper [41] by making use of a scalar-tensor theory of gravity and the 
principles of the analytic theory of relativistic reference frames in the solar system [63, 75, 76] that 
was adopted by the XXIV-th General Assembly of the International Astronomical Union [5, 77] as 
a standard for data processing of high-precision astronomical observations. 

The advanced post-Newtonian dynamics of the Sun-Earth-Moon system must include the fol- 
lowing structural elements: 

1. construction of a set of astronomical reference frames decoupling orbital dynamics of the 
Earth-Moon system from the rotational motion of the Earth and Moon with the full account 
of the post-Newtonian corrections and elimination of the gauge modes; 

2. relativistic definition of the integral parameters like mass, the center of mass, the multipole 
moments of the gravitating bodies; 

3. derivation of the relativistic equations of motion of the center-of-mass of the Earth-Moon 
system with respect to the barycentric reference frame of the solar system; 
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4. derivation of the relativistic equations of motion of the Earth and Moon with respect to the 
reference frame of the Earth-Moon system; 

5. derivation of the relativistic equations of motion of CCR on the Moon (or a lunar orbiter 
that is deployed with CCR) with respect to the selenocentric reference frame; 

6. derivation of the relativistic equations of motion of a laser with respect to the geocentric 
reference frame. 

These equations must be incorporated to LLR data processing software operating with observable 
quantities, which are proper times of the round trip of the laser pulses between the laser on the 
Earth and CCR on the Moon. The computational advantage of the new approach to the lunar 
ephemeris is that it separates clearly physical effects from the choice of coordinates. This allows 
us to get robust measurement of true physical parameters of the LLR model and give them direct 
physical interpretation. The new approach is particularly useful for comparing different models of 
the lunar interior and for making the fundamental test of general theory of relativity. 
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